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ABSTRACT 


This  paper  is  concerned  with  the  propagation  of  phase  boundaries 
in  elastic  bars.  It  is  known  that  the  Riemann  problem  for  an  elastic 
bar  capable  of  undergoing  isothermal  phase  transitions  need  not 
have  a  unique  solution,  even  in  the  presence  of  the  requirement 
that  the  entropy  of  any  particle  cannot  decrease  upon  crossing  a 
phase  boundary.  For  a  special  class  of  elastic  materials,  we  have 
shown  elsewhere  that,  if  all  phase  boundaries  move  subsonically 
with  respect  to  both  phases,  this  lack  of  uniqueness  can  be  resolved 
by  imposing  a  nucleation  criterion  and  a  kinetic  relation  for  the 
relevant  phase  transition.  Here  we  investigate  an  alternative 
approach  that  singles  out  acceptable  solutions  on  the  basis  of  a 
theory  that  adds  effects  due  to  viscosity  and  second  strain  gradient 
to  the  elastic  part  of  the  stress.  We  show  that,  for  phase  boundaries 
that  propagate  subsonically,  this  approach  is  equivalent  to  the 
imposition  of  a  particular  kinetic  relation  at  the  interface  between 
the  phases. 
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bar  according  to  a  one-dimensional  continuum  model  that  permits  the  material  of  the  bar  to 
undergo  phase  transitions.  Although  inertial  effects  were  taken  into  account  in  [1],  motions 
were  assumed  to  take  place  isothermally.  When  a  motion  of  the  bar  involves  a  propagating  strain 


discontinuity,  the  second  law  of  thermodynamics  imposes  an  entropy  admissibility  requirement: 
at  each  strain  jump,  the  product  of  the  velocity  of  the  discontinuity  and  an  associated  driving 
traction  must  be  non-negative. 

If  the  material  of  the  bar  is  such  that  stress  is  a  montonically  increasing,  strictly  convex  or 
strictly  concave  function  of  strain,  then  phase  transformations  cannot  occur,  and  all  propagating 
strain  discontinuities  are  shock  waves.  For  a  bar  made  of  such  a  material,  the  Cauchy  problem 
for  the  associated  field  equations  and  jump  conditions  has  at  most  one  solution  that  fulfills  the 
entropy  admissibility  requirement. 

On  the  other  hand,  for  the  class  of  materials  considered  .r  [1],  the  stress-strain  relation  is 
such  that  stress  at  first  increases  with  increasing  strain,  then  decreases  and  finally  increases 
again.  A  constitutive  law  of  this  kind  affords  a  continuum  model  for  phase  transitions  in  the 
present  setting.  A  propagating  strain  discontinuity  may  now  be  either  a  shock  wave  or  a  phase 
boundary,  according  to  whether  the  particles  separated  by  the  discontinuity  are  of  the  same  phase 
or  of  distinct  phases.  In  this  setting,  the  Cauchy  problem  need  no  longer  have  a  unique  solution, 
even  with  the  entropy  admissibility  requirement  in  force.  For  the  Riemann  problem,  it  was 
shown  in  [1]  that  this  lack  of  uniqueness  may  be  viewed  as  arising  from  the  need  to  specify  two 
additional  pieces  of  constitutive  information  pertaining  to  phase  boundaries:  a  nucleation 
criterion  for  the  initiation  of  a  phase  transition  and  a  kinetic  relation  that  controls  the  rate  at 
which  the  phase  transition  proceeds.  The  importance  of  nucleation  and  kinetics  has  long  been 
recognized  in  the  materials  science  literature  concerning  phase  transformations  [2]. 

A  similar  need  had  been  identified  in  an  earlier  study  of  quasi-static  processes  involving 
phase  transformations  in  elastic  bars  [3].  The  model  for  such  processes  given  in  [3]  incorporates 
a  nucleation  criterion  and  a  kinetic  relation.  It  leads  to  a  determinate  quasi-static  theory  and 
yields  results  that  are  in  qualitative  agreement  with  experimental  observations  for  bars  made  of 
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"shape-memory"  materials  [4].  The  form  of  the  kinetic  relation  employed  in  [3]  is  one  in  which 
the  driving  traction  is  a  materially-determined  function  of  the  velocity  of  the  moving  phase 
boundary;  the  nucleation  criterion  involves  the  specification  of  a  critical  level  of  driving  traction 
that  signals  the  onset  of  a  phase  transformation. 

For  a  special  material  whose  rising-falling-nsing  stress-strain  curve  is  piecewise  linear, 
we  showed  in  [1]  that,  for  the  Riemann  problem,  the  extent  of  lack  of  uniqueness  of  solution 
remaining  after  imposition  of  the  entropy  admissibility  requirement  was  precisely  that  needed  to 
accommodate  the  nucleation  criterion  and  the  kinetic  relation  at  phase  boundaries,  at  least  for 
kinetic  relations  of  the  kind  introduced  in  [3]  and  in  circumstances  where  the  phase  boundary 
propagates  subsonically  with  respect  to  both  phases.  Kinetic  relations  may  not  be  prescribed  at 
shock  waves.  Phase  boundaries  that  move  supersonically  with  respect  to  the  material  phase  with 
lower  sound  speed  were  also  discussed  in  [1];  it  was  found  that  no  kinetic  relation  may  be 
prescribed  at  such  a  phase  boundary  in  the  Riemann  problem. 

A  different  approach  to  the  issue  of  identifying  meaningful  solutions  to  a  quasilinear 
system  of  dynamical  field  equations  and  jump  conditions  in  the  presence  of  propagating  phase 
boundaries  has  been  put  forward  by  Slemrod  [5,6,7]  and  explored  by  Shearer  [8,9,10];  related 
ideas  have  been  pursued  by  Truskinovsky  [11,12]  and  others.  Although  these  authors  are  often 
concerned  with  the  van  der  Waals  fluid,  they  note  that  their  method  can  be  applied  to  elastic  bars 
as  well.  In  this  approach,  one  first  augments  the  conventional  constitutive  law  for  an  elastic 
material  capable  of  undergoing  isothermal  phase  transitions  in  such  a  way  that  the  stress  a 
depends  not  only  on  the  strain  y,  but  also  on  the  strain  rate  y{  and  the  second  spatial  strain 
gradient  y  .  The  dependence  of  a  on  y  is  nonlinear,  but  y  and  y  enter  linearly  and  are 
associated  with  a  viscosity  v  and  a  strain-gradient  coefficient  X,  respectively.  Next,  one 
establishes  a  criterion  to  identify  those  moving  phase  boundaries  within  the  elastic  theory  that 
are  preferred  on  the  basis  of  this  augmentation:  the  criterion  asserts  that  such  a  phase  boundary  is 
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preferred  if  and  only  if  the  strains  y  and  on  either  side  of  the  discontinuity  can  be  smoothly 
connected  by  a  traveling  wave  constructed  within  the  augmented  theory. 

The  present  paper  is  concerned  with  the  relationship  between  the  approach  based  on  such 
an  augmented  model  and  the  direct  approach  described  in  [1]. 

In  the  next  section,  we  formulae  the  elastic  theory  and  its  augmentation  as  described 
above.  In  Section  3,  we  review  the  local  properties  of  strain  discontinuities  according  to  the 
elastic  theory;  as  in  [1],  we  impose  a  kinetic  relation  at  phase  boundaries.  We  consider  the 
augmented  theory  in  Sections  4-6,  where  we  study  traveling  waves  in  the  manner  of  Slemrod 
[5,6,7],  but  for  the  case  in  which  the  strain  dependence  of  stress  is  trilinear.  We  show  in  Section 
5  that  strains  y  and  y  can  be  connected  by  a  subsonic  traveling  wave  in  the  augmented  theory  if 
and  only  if  they  are  related  by  a  certain  restriction.  In  contrast,  no  such  restriction  applies  to 
strains  that  are  connectable  by  a  supersonic  traveling  wave;  see  Section  6.  In  Section  7,  we 
delineate  the  class  of  phase  boundaries  within  the  elastic  theory  that  are  preferred  by  the  criterion 
described  above.  In  the  subsonic  case,  we  find  that  the  preferred  phase  boundaries  are  those  that 
conform  to  a  particular  kinetic  relation  that  arises  directly  from  the  restriction  mentioned  above. 
This  kinetic  relation  is  of  precisely  the  kind  introduced  in  the  quasi-stat'c  setting  in  [3]  and 
studied  in  the  dynamic  context  in  [1],  For  supersonic  phase  boundaries,  no  kinetic  relation 
emerges  from  the  augmented  theory.  The  distinction  between  subsonic  and  supersonic  phase 
boundaries  observed  here  is  entirely  consistent  with  the  findings  reported  in  [1]  and  described 
briefly  above. 

2.  Basic  equations.  Suppose  a  bar  of  unit  cross-sectional  area  occupies  the  interval 
(-oo,  oo)  of  the  x-axis  in  a  reference  configuration.  During  a  motion,  the  particle  at  x  in  the 
reference  configuration  is  carried  to  x  +  u(x,t)  at  time  t,  where  u(x,t)  is  the  diplacement.  We 
assume  that  u  is  at  least  continuous  with  piecewise  continuous  first  and  second  derivatives 
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throughout  the  regions  of  space-time  to  be  considered;  more  stringent  requirements  will  often  be 
needed.  The  strain  y(x,t)  and  particle  velocity  v(x,t)  are  defined  by 


y=ux  ,  v  =  ut,  (2.1) 

wherever  the  derivatives  exist.  Balance  of  momentum  and  compatibility  of  y,v  require  that 
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at  points  where  cr,  v,  and  y  are  smooth;  here,  c(x,t)  is  the  stress  at  (x,t),  and  p  is  the  (constant) 
mass  density  in  the  reference  configuration.  If  either  y  or  v  is  discontinuous  across  a  curve 
x  =  s(t)  in  the  -plane,  balance  of  momentum  and  the  smoothness  properties  of  u  yield  the 
following  jump  conditions: 


[[<?]]  =  -ps  [[v]]  (2.4) 

[[v]]  =  -s[[y]],  (2.5) 

where  s(t)  is  the  velocity  of  the  moving  strain  discontinuity,  and  for  any  g(x,t)  we  have  written 
[[g]]  =  g(s(t)+,t)  -  g(s(t)-,t). 

For  the  elastic  bar,  we  take 


—  _  A-/.  A 
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here  &(y)  is  the  stress  response  function  whose  graph  is  shown  in  Figure  1.  The  strain  y  is 
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restricted  to  the  range  (-1,  +°°)  in  order  to  assure  that  the  mapping  x  ->  x  +  u(x,t)  is  invertible  for 
each  t.  We  say  that  a  particle  of  the  bar  labeled  by  x  is  in  phase  1,2  or  3  at  time  t  according  to 
whether  the  strain  y(x,'0  is  in  (-1,  yM],  (yM,  ym)  or  [ym,  <*>),  respectively.  The  constants  |i,j  and  p. 

<  (ij  are  the  elastic  moduli,  and  Cj  =  (fij/p)  and  c3  =  (p3/p)  <  Cj  are  the  corresponding 

sound  wave  speeds,  in  phases  1  and  3.  The  special  stress  gq  shown  in  Figure  1  is  such  that  the 

two  shaded  areas  are  equal;  it  is  called  the  Maxwell  stress.  For  the  trilinear  material, 

1/2 

a  =  (a.  ,o  )  ,  where  a.,  and  G  are  the  respective  stresses  at  the  local  maximum  and 
o  '  M  nr  M  in  r 

A 

minimum.  On  the  declining  branch  of  the  stress-strain  curve,  a(y)  -  -  p2y  +  c^.  We  set 

1/2  a 

c2=(p2/p)  ;  note  that  c2  is  not  a  wave  speed.  Continuity  of  o(y)  at  the  maximum  and  the 

minimum  requires  that 


o2/P  =  (c^4)ym  =  (c^c^ 


m 


(2.7) 


In  the  augmented  theory,  the  constitutive  statement  (2.6)  is  replaced  by 


G  =  o(y)  +  pvyt-pXyxx  , 


(2.8) 


where  G(y)  is  again  the  trilinear  function  of  Figure  1,  and  v  >  0  and  X  >  0  are  the  viscosity  and 
the  strain-gradient  coefficient,  respectively.  The  dimensionless  parameter  co  def  ned  by 


co  =  2  Xl,2/\ 


(2.9) 


will  be  useful  in  what  Allows. 

When  considering  the  material  governed  by  (2.8),  we  shall  str  ngthen  the  smoothness 
requirements:  i.  will  be  twice  continuously  differentiable  with  piecewise  continuous  third  and 
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fourth  derivatives  in  all  parts  of  space-time  to  be  considered.  Of  the  jump  conditions  (2.4),  (2.5), 
the  latter  is  then  trivially  satisfied,  and  the  former  reduces  to  the  requirement  that  a  be 
continuous;  by  (2.8),  this  requires  that  y  ,  and  therefore  u  ,  be  continuous. 

XX  XXX 

The  augmented  theory  used  here  is  analogous  to  a  special  case  of  a  theory  proposed  by 
Korteweg  [13]  to  account  for  the  effect  of  capillarity  in  a  viscous  fluid;  see  [5].  In  the  latter 
theory,  one  identifies  y  in  (2.8)  with  specific  volume,  -a  with  pressure,  p  with  density,  v  with 
viscosity  and  A.  with  capillarity.  Thermo-mechanical  difficulties  raised  by  Korteweg’ s  theory 
have  been  discussed  in  [14]. 

3.  Local  theory  of  phase  boundaries  in  elastic  bars.  We  now  consider  the  elastic  case 
(2.1)  -(2.6),  together  with  the  appropriate  smoothness  assumptions.  Although  many  of  the 
detailed  results  to  follow  apply  only  for  the  triiinear  material,  the  more  general  discussion  in  the 
first  part  of  the  present  section,  as  well  as  certain  subsequent  results,  are  valid  for  any 
rising-falling-rising  stress-strain  curve. 

Suppose  that,  during  a  motion,  the  portion  of  the  bar  occupying  the  interval  [Xj,  x2]  in  the 
reference  configuration  bears  fields  y,  v  and  o  that  are  smooth  except  for  a  single  moving  strain 
jump  at  x  =  s(t).  Let 


a(y')dy' 


y>-l 


0 


(3.1) 


be  the  stored  energy  per  unit  volume  of  the  bar,  and  consider  the  total  mechanical  energy 
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(3.2) 


1 

W(y)  +  -  pv 


1 


dx 


x 


1 


associated  with  the  piece  of  the  bar  under  consideration  at  time  t.  A  direct  calculation  using 
(3.1),  (2.2) -(2.5)  shows  that 

c(x2,t)v(x2,t)  -  o(xltt)v(xltt)  -  E(t)  =  f(t)'s(t) ,  (3.3) 


where  the  driving  traction  f  acting  on  the  strain  discontinuity  is  defined  by 


f=f(Y,Y)=J  o(Y)dY  -  ^(i  +  cj)(f -y) 


(3.4) 


Y 


+  +  +  + 

Here  Y  =  Y(0  =Y(s(t)±,t)  and  a  =  a(t)  =  a(s(t)±,t)  are  the  strains  and  stresses,  respectively,  on  the 
positive  and  negative  sides  of  the  strain  discontinuity.  The  left  side  of  (3.3)  represents  the  excess 
of  the  rate  of  work  of  the  external  forces  over  the  rate  of  increase  of  mechanical  energy. 

Although  this  excess  vanishes  for  smooth  strain  fields  (f=0)  or  for  stationary  discontinuities 
(s=0),  this  is  not  the  case  in  general.  If  the  material  is  viewed  as  thermoelastic,  and  if  we  make 
the  assumption  -  however  unrealistic  -  that  the  motions  considered  take  place  isothermally,  then 
as  shown  in  [15],  the  Clausius-Duhem  inequality  requires  that  the  instantaneous  rate  of 
mechanical  dissipation  f(t)s(t)  be  non-negative: 

f(t)s(t)  >  0  (3.5) 

for  all  time.  The  entropy  admissibility  condition  (3.5)  must  hold  at  all  strain  discontinuities; 
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under  isothermal  conditions,  it  is  equivalent  to  the  assertion  that  the  entropy  of  a  particle  cannot 
decrease  as  the  particle  crosses  a  discontinuity  in  strain. 

At  the  moving  discontinuity  x=s(t),  the  jump  conditions  (2.4),  (2.5)  imply 


A  4.  A  . 

o(y)  -  g(y) 


?-Y 


(3.6) 


[a(y)  -  o(y)](y  -  Y)  =  P  (v  -  v)2  . 


(3.7) 


Tne  right  side  of  (3.6)  is  thus  necessarily  non-negative  for  any  pair  of  strains  Y>  Y  that  can  occur 
at  a  strain  jump. 

Conversely,  if  y,  Y  are  numbers  in  (-1,  °°)  such  that  the  right  side  of  (3.6)  is  non-negative, 

-  +  • 

then  it  is  possible  to  find  numbers  v,  v  and  s  such  that  the  jump  conditions  (2.4),  (2.5)  are 
satisfied. 


Because  we  are  concerned  here  with  phase  boundaries  rather  than  shock  waves,  we 
assume  that  y  and  Y  belong  to  different  phases.  Since  we  shall  not  consider  here  cases  in  which 

-  +  4. 

either  y  or  y  is  in  phase  2,  i.e.,  the  "unstable  phase",  we  may  as  well  take  y  to  be  in  phase  1,  and  y 
in  phase  3.  When  specialized  to  this  case  and  to  the  trilinear  material  of  Figure  1,  (3.6)  becomes 


2  +  2- 
•2  =  SY^i_Y 

f-Y 


(3.8) 


I  -4-  -  4- 

In  the  y,  Y-plane,  the  set  ot  pairs  y,  y  for  which  y  is  in  phase  1,  y  is  in  phase  3  and  the 
right  side  of  (3.8)  is  non-negative  is  represented  by  the  shaded  region  T  shown  in  Figure  2.  At 
any  point  on  the  boundary  segment  BC,  the  numerator  in  (3.8)  vanishes,  so  that  s  =  0.  The 


I 

corresponding  phase  boundary  is  thus  instantaneously  stationary.  In  the  special  case  of  | 

time-independent  fields,  points  (y,  *y0  on  BC  correspond  to  equilibrium  phase  mi  ;tures.  If  oQ  j 

denotes  the  equal-area  stress  in  Figure  1,  the  special  equilbrium  mixture  for  which  y  =  oJpj  ,  | 

g 

+  .  3 

y  =  Gy/)i3  is  called  the  Maxwell  state',  it  corresponds  to  point  M  in  Figure  2.  'Hie  Maxwell  state  jj 

* 

represents  stable  equilibrium;  other  phase- 1,  phase-3  equilibria  are  merely  metastable.  \ 


The  driving  traction  f  acting  on  a  phase- 1,  phase-3  interface  can  be  found  from  (3.4)  and 

A 

the  explicit  form  of  o(y)  for  the  trilinear  material:  one  finds  that 


It  follows  that  the  driving  traction  vanishes  on  the  hyperbola  y  f  =  yMym  .  In  view  of  (3.3),  this 
means  that  points  in  T  that  lie  on  this  hyperbola  correspond  to  values  of  y  and  y  for  which  the 
associated  phase  boundary  x  =  s(t)  propagates  without  dissipation.  Note  that  the  Maxwell  point 
M  lies  on  the  hyperbola.  At  points  off  the  hyperbola,  f  *  0,  so  that  the  entropy  admissibility 
requirement  (3.5)  determines  the  sign  of  s;  see  Figure  2.  For  points  on  the  hyperbola  f  =  0,  the 
sign  of  s  is  not  determined  by  (3.5),  and  propagation  in  either  direction  is  possible. 

With  the  help  of  (3.8)  and  (3.5),  one  can  show  (as  in  [1])  that  s  necessarily  lies  in  the 

range 


-c3  <  s  <  c*  , 


(3.10) 


where 


-10- 


(3.11) 


c*  = 


V3 


,1/2 


1  +7 


C3< 


C*< 


m 


The  propagation  speed  s  is  said  to  be  subsonic  if  is!  <  c3  .  From  (3,8),  it  follows  that  this  occurs 
only  for  those  points  in  the  region  T  of  Figure  2  for  which  y  >  0.  For  supersonic  motion  of  the 
phase  boundary,  y  must  be  negative,  so  that  part  of  the  bar  is  in  compression.  Furthermore,  the 
driving  traction  is  positive,  so  that  the  entropy  admissibility  requirement  forces  the  direction  of 
propagation  of  the  phase  boundary  to  be  such  that  the  portion  of  the  bar  in  phase  1  grows  at  the 
expense  of  phase  3.  Thus  for  a  supersonically  propagating  phase  boundary,  the  "parent"  phase  - 
the  phase  into  which  the  interface  moves  -  is  always  phase  3. 

We  now  use  (3.8)  and  (3.9)  to  map  the  region  T  of  the  y,  '/-plane  into  the  s, f-plane.  Each 
point  (y,  $)  that  does  not  lie  on  BC  (Figure  2)  is  carried  to  two  points  (s,  f)  and  (-s,  f)  in  the 
s,  f-plane;  if  f  *  0,  only  one  of  these  satisfies  the  entropy  admissibility  requirement  (3.5).  If  f=  0, 
the  point  (y,  y)  lies  on  the  hyperbola;  such  a  point  maps  to  a  pair  of  admissible  points  (±s,  0)  in 
the  s,  f-plane,  with  0  <  s  <  c3  .  Each  point  on  BC  maps  to  a  single  admissible  point  (0,  f)  on  the 
f-axis.  All  points  in  T  that  lie  on  y=0  map  to  the  same  pair  of  points  (±c3,  fQ),  where  the  constant 
fQ  is  given  by  fQ=(|J-1  -  M-3)  YmYM/2  >  0.  Figure  3  shows  the  admissible  image  of  T  in  the 
s,  f-plane. 

The  physical  basis  for  the  elucidation  of  phase  transitions  in  solids  involves  both  a 
nucleation  criterion  governing  the  initiation  of  the  transition  and  a  kinetic  relation  controlling  the 
rate  at  which  it  proceeds;  see  [2, 16, 17].  We  have  discussed  simple  continuum-mechanical 
implementations  of  these  notions  in  [1,  3, 15].  In  particular,  it  is  shown  in  [1]  that,  as  long  as  all 
phase  boundaries  propagate  subsonically,  a  nucleation  criterion  and  a  kinetic  relation  can  be 
imposed  in  the  Riemann  problem,  and  that  they  single  out  unique  admissible  solutions  in  that 
context.  In  the  present  discussion,  we  are  concerned  only  with  the  propagation  of  an  existing 
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phase  boundary  and  not  w;,h  its  emergence;  the  nucleation  criterion  is  therefore  not  relevant 
here.  As  for  the  kinetic  relation,  we  proceed  as  in  [  1]  and  assume  that,  if  the  phase  boundary 
velocity  s  is  subsonic,  there  is  a  function  <p  determined  by  the  material  that  relates  the  driving 
traction  f  to  s,  the  latter  being  a  measure  of  the  rate  at  which  the  phase  transition  takes  place. 
Thus  we  take 


f  —  cp(s)  ,  -c3<s<c3  .  (3.12) 

Because  of  the  entropy  admissibility  requirement  (3.5),  cp  must  satisfy 

cp(s)s > 0  ,  -c3<s<c3.  (3.13) 

Note  that  if  tp(s)  is  continuous  at  s  =  0  as  we  shall  assume  here,  (3.13)  requires  that  tp(0)  =  0.  It  is 
required  that  the  curve  represented  by  (3.12)  lie  in  the  hatched  region  of  the  s,f-plane  shown  in 
Figure  3.  The  pre-image  of  this  curve  in  the  y,  y-plane  under  the  mapping  (3.8),  (3.9),  (3.5)  thus 
comprises  the  locus  of  all  strain-pairs  y,  f  at  a  phase- 1,  phase-3  interface  that  are  consistent  with 
the  kinetic  relation  (3.12). 

The  choice  cp(s)  =  0  for  the  kinetic  response  function  corresponds  to  dissipation-free 
motions  of  the  bar.  For  nearly  dissipation-free  motions,  the  choice 

<p(s)  =  Ks  ,  -c3<s<c3,  (3.14) 


..  ~  ~ , rr ! - : *.i. .  , _ n  _ „ 

Vvliu  is.  a  bUiucAcnuy  Miicui  liOn-ncgauvc  cuiiduuu,  io  iiaiuicu. 


- 1 /o  i  a\ 

l'wiiiV/Ui’  luiatiuiio  uuiw  man 


can  of  course  be  considered.  The  choice 
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<p(s)  =  Ksinh’^s/k)  , 


(3.15) 


where  K  and  k  are  constants,  can  be  motivated  by  arguments  of  the  kind  used  in  thermal 
activation  theory;  for  a  sketch  of  the  latter  theory,  see  Chapter  1  of  [16].  Other  special  choices  of 
cp(s)  are  discussed  in  [1,3].  Kinetic  relations  not  of  the  form  (3.12)  may  be  of  interest.  One  that  is 
suggested  by  a  dislocation-based  theory  put  forward  in  [17]  for  the  modeling  of  fast  martensitic 
phase  transitions  is 


f  =  f„  +  ms  +  ks  , 


(3.16) 


with  f*  ,  m  and  k  constants. 

The  kinetic  relations  discussed  above  represent  local  requirements,  as  indeed  they  must  if 
they  are  to  be  regarded  as  part  of  the  constitutive  description  of  the  material.  Moreover,  the 
notion  of  a  scalar  driving  traction,  and  consequently  that  of  a  kinetic  relation,  can  be  carried  over 
to  more  general  constitutive  settings  in  three  dimensions  with  thermal  effects  accounted  for  as 
well;  for  further  discussion,  see  [15]. 

It  was  shown  in  [1]  that  the  Riemann  problem  based  on  the  theory  of  elastic  bars 
described  above  is  not  overdetermined  when  the  kinetic  relation  (3.12)  is  imposed  at 
subsonically  propagating  phase  boundaries. 

4.  The  augmented  theory:  construction  of  solutions.  In  order  to  pursue  the  approach 
based  on  the  augmented  theory,  we  replace  the  constitutive  law  (2.6)  by  (2.8),  and  we  replace 

4*  4* 

the  moving  discontinuity  with  strains  y  and  particle  velocities  v  in  the  elastic  bar  by  a  traveling 
wave  in  the  augmented  theory.  In  the  traveling  wave, 
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7  =  7(4) »  v  =  v®,  $  =  x-st, 


(4.1) 


where  s  is  a  constant,  and 


7(-°°)  =  7  *  v(-°°)  =  v ,  7(+<»)  =  $  ,  v(+oo)  =  v  ,  (4.2) 

corresponding  to  a  smooth  connection  of  the  states  7,  v  and  $ ,  v.  As  before,  7  is  to  be  in  phase  1, 
^  in  phase  3.  In  conformity  with  the  smoothness  assumptions  made  in  Section  2,  we  suppose 
that  y(t;)  and  v(%)  are  continuous  and  continuously  differentiable  on  (-<*>,  °°),  and  that  they  have 
piecewise  continuous  second  and  third  derivatives  there.  In  view  of  (2.4)  and  (2.8),  y  "(£)  must 
in  fact  be  continuous  for  all  £,,  provided  kt 0. 

From  (2.2),  (2.3)  and  (2.8),  traveling  waves  of  the  form  (4.1)  must  satisfy 

pXy  +  pvs  7"-  a'(y)  7'  -  ps  v '  =  0 ,  (4.3) 

v'  +  sy'  =  0,  (4.4) 

at  all  points  where  y  '"(£)  exists.  We  seek  solutions  7,  v  of  (4.3),  (4.4)  that  satisfy  the  conditions 
(4.2)  at  ±o«. 

Before  constructing  an  explicit  solution  to  this  boundary  value  problem,  we  record  some 
useful  observations.  Assuming  the  existence  of  a  solution,  we  first  integrate  (4.3),  (4.4)  to  obtain 

.  „  .,1a 

Xy  +vsy  --c?(y)-sv  =  A,  (4.5) 

P 

v  +  sy=B,  (4.6) 

where  A  and  B  are  constants.  From  the  boundary  conditions  (4.2),  it  follows  that  A  and  B  satisfy 
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(4.7) 


1  A  +  .4.  1  A  .  .. 

A  =  -  -  o(Y)  -  s  v  =  -  -  c(y)  -  s  v  , 
P  P 


B  =  v  +  sf  =  v  +  sy. 


But  (4.7)  immediately  imply  (3.6)  and  (3.7),  which  are  therefore  necessary  conditions  on  s  and 
the  given  data  for  the  existence  of  a  solution  to  the  problem  (4.2)-(4.4). 

Elimination  of  v  between  (4.3)  and  (4.4)  yields  a  differential  equation  for  y: 

X  y  "  +  vs  y '  -  -  o(y)  +  s2  y  =  C  on  (-oo,  <*>) ,  (4.8) 

P 


where  the  constant  C  is  given  by 


.  <!t  .i  +  a  .  o . 

C  =  A  +  sB  =  -~+s  y  =  --+s  Y- 
P  P 


(4.9) 


41  ^  4" 

and  we  have  written  a  s  o(y).  The  boundary  conditions  associated  with  (4.8)  are 


Y(- °°)  =  Y.  y(+~)  =  f,  (4.10) 

with  yin  phase  1,  y  in  phase  3.  We  speak  of  (4.8)-(4.10)  as  the  traveling  wave  problem.  Once 
y(£,)  has  been  found,  v(£,)  is  determined  from  (4.6). 

After  multiplying  the  differential  equation  (4.8)  by  y  '(£)  and  integrating  with  respect  to  E, 
from  -  oo  to  +  oo,  one  finds  with  the  help  of  (3.6),  (4.9)  that  a  solution  y(£),  if  it  exists,  must 
satisfy 
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oo 


(4.11) 


f(Y>'?)  =  pvs  J  [YU)]2d^, 

-oo 


where  f  is  the  function  defined  in  (3.4).  From  (4.1 1)  and  the  fact  that  v  >  0,  it  follows  that,  if  the 
traveling  wave  problem  has  a  solution,  then  s  and  the  data  y,  y  necessarily  satisfy 

f  (y,  y)  s  >  0 ;  (4.12) 

cf.  the  entropy  admissibility  condition  (3.5)  of  the  elastic  theory.  Of  the  two  necessary  conditions 

(3.6)  and  (4.12),  the  first  determines  the  speed  of  propagation  Isl  of  the  traveling  wave.  The  sign 

of  s,  and  hence  the  direction  of  propagation,  is  then  determined  by  (4.12)  with  the  help  of  (3.4), 

±  A  -  + 

provided  the  given  strains  y  are  such  that  f (y,  y)  *  0.  As  in  the  derivation  of  (3.10),  from  (3.6) 
and  (4.12)  one  finds  that  s  lies  in  the  interval  (-c3,  c„). 

We  seek  a  solution  of  the  traveling  wave  problem  (4.8)-(4.10)  for  which  y(£)  is  in  phase 
1  for  -oo  <  ^  <  0,  in  phase  2  for  0  <  £  <  b,  and  in  phase  3  for  b  <  2;  <  °° ,  for  some  suitable  b  to  be 
determined.  In  view  of  our  smoothness  assumptions,  the  following  conditions  must  hold  at  the 
interfaces  between  the  various  phases: 

Y(°-)  =  Ym  >  Y(0+)  =  Ym  .  Y  '(0-)  =  Y  '(0+) ,  (4.13) 

Y(t>-)  =  Ym  .  Y(b+)  =Ym  >  Y  '(•>)  =  Y  '(b+) .  (4.14) 

If  (4.13),  (4.14)  and  the  differential  equation  (4.8)  all  hold,  then  y"(£)  is  automatically 
continuous  at  H.-Q  and  %=b,  provided  X*Q. 

In  order  to  construct  such  a  solution,  we  first  consider  three  subsidiary  problems  for  the 
restrictions  of  Yfo  the  separate  intervals  (-<*>,  0],  [0,  b]  and  [b,  °°).  The  first  of  these  is  a 
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boundary-value  problem,  while  the  second  and  third  are  initial-value  problems;  the  problems  are 
to  be  solved  in  order. 

Problem  1: 

Xy"  +  vsy/-(Cj -s2)y=C  on 
Y(-°°)  =  7 .  7(0-)  ==7;^  ; 

Problem  2: 

X  7  "  +  vs  7 '  +  (c2  +  s2)  7  =  C  +  a2/p  on  [0,  b] , 

7(°+)  =  7m*  7  '(0+)  =  7  '(0-)  ; 

Problem  3: 

A,  7  "  +  vs  7 '  -  (c2  -  s2)  7  =  C  on  [b,  °°) , 

7(b+)=7m.  7'(b+)  =  7'(b-) . 

The  function  y  on  (-<»,  determined  by  solving  these  problems  will  satisfy  the  differential 
equation  (4.8),  the  first  of  the  boundary  conditions  (4.10),  all  of  the  interface  conditions  (4.13)  at 
£=0,  and  all  but  the  first  of  the  interface  conditions  (4, 14)  at  £=b.  Thus  after  solving  Problems 
1-3,  we  shall  still  need  to  enforce  the  boundary  condition 

Y(+oo)  =  y  ,  (4.21) 


(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 

(4.20) 


the  interface  condition 


Y(b-)  =  Ym 


(4.22) 
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and  the  requirements 


-1  <Y(5)^YM  for-°c<^<0;  YM  <7(4)  <  Ym  for  0  <  £,  <  b;  y®  >  ym  for  E,  >  b .  (4.23) 

The  inequalities  (4.23)  insure  that  the  strain  belongs  to  the  appropriate  phase  in  each  interval. 

After  using  the  right-most  representation  of  C  in  (4.9),  one  finds  that  the  solution  of 
Problem  1  is  given  by 

Y(^)=Y+(YM-Y)ePl^  ,  -°o<£<0,  (4.24) 

where  the  positive  number  p}  is  defined  by 

p!  ^  (s  - p)  >  °  ,  P  =  P(s)  =  [s2  +  co2  (c2  -  s2)]1/2  >  0 ,  (4.25) 

and  co  is  defined  in  (2.9). 

We  turn  next  to  Problem  2.  The  general  solution  of  (4.17)  is 

7(5)  =  Y0  +  D1eq^  +  D2 ,  0  <  %  <  b  ,  (4.26) 

where  D1  and  D2  are  arbitrary  constants,  and  the  constant  is  defined  by 

Y0  =  (C  +  a2/p)/(c2  +  s2)  .  (4.27) 

In  (4.26),  the  (possibly  complex)  constants  q  and  q2  in  the  exponentials  are  given  by 
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ql  '  2X  S "  Q  ’  q2  2X (S  +  Q)  ’ 


(4.28) 


where  we  have  set 


Q  =  Q(s)  = 


r-2  2/2  •  2m  1/2  -r  -2  2/2  -2^/, 

[s  -0)  (c2  +  S  ))  '  if  S  -CO  (c2  +  s  )>0, 

.  r  2/2  -2X  - 2t  1/2  -2  2/2 

1  [CO  (c2  +  S  )  -  S  ]  if  S  -  CO  (c2  +  S  )  <  0  . 


(4.29) 


We  now  choose  the  constants  Dj  and  D2  so  that  the  initial  conditions  (4.18)  for  Problem  2  are 
satisfied.  This  leads  to 


D  =  - 


(Y, 


M  Vp^P2' V  ’  °2  %Q^M-Vp2(p2-q2>; 


(4.30) 


here  p2  is  given  by 


p2  =  -^(s  +  p)<°, 


(4.31) 


and  P  is  defined  in  (4.25).  It  has  been  assumed  above  that  Q*0,  so  that  »  while  the  case 
Q=0  can  be  readily  treated,  we  exclude  it  for  simplicity.  Note  that  although  and  q2  may  be 
either  real  or  complex,  y  as  given  on  [0,  b]  by  (4.26)-(4.31)  is  always  real. 


Finally,  we  consider  Problem  3.  Its  solution  is  given  by 


1© 


+  V  f  r  r,£  -  b>  .  r9^  -  b)l  .  Y'(b-)  r  r,(g  -  b)  r,($  -  b)-]  c 

=  y  +  — —  -  r2  e  1  +  Tje  L  + -  e  1  ,  q  >  b ,  (4.32) 

rl  -  r2  L  J  Tj  -  r2  L  J 
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where  r ^  and  r2  are  defined  by 


r 


l 


(4.33) 


and  we  have  set 


R  =  R(s)  = 


[s2  +  co2(C2-k2)]1/2  if  s2  +  co2  (c2  -  s2)  >  0  , 
i  [co2  (s2  -  c2)  -  s2]1/2  if  s2  +  co2  (c2  -  s2)  <  0 . 


(4.34) 


In  (4.32),  y  '(b-)  is  found  from  (4.26): 

y  '(b-)  =  Djqj  eqib  +  D2q2  eQ2b  ,  (4.35) 

with  Dj  and  D2  given  by  (4.30).  For  simplicity,  we  have  excluded  the  case  R=0,  corresponding 
to  Tj=  r2-  Again,  we  note  that,  although  r}  and  r2  may  be  either  real  or  complex,  y(^)  as  given  on 
[b,oo)  by  (4.32)-(4.35),  (4.28)-(4.30)  is  always  real. 

In  order  to  complete  the  construction  of  a  solution  y (^)  to  the  traveling  wave  problem 
from  the  solutions  of  Problems  1, 2  and  3,  it  is  necessary  to  address  the  remaining  requirements 
(4.2'i)-(4.23).  By  (4.26),  we  may  write  (4.22)  as 

D,  eqlb  +  D„  eq2b  =  y  -yn,  (4.36) 

i  z.  iTi  \j 

where  Dj  and  D2  are  given  by  (4.30),  q}  and  q2  by  (4.28)  and  yQ  by  (4.27).  Equation  (4.36)  is  a 
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condition  for  the  determination  of  the  length  b  of  the  phase-2  interval;  we  postpone  the 
discussion  of  the  solvability  of  (4.36). 

Of  the  requirements  (4.23),  the  first  pertains  to  the  interval  (-«»,  0];  by  (4.24),  this 
condition  is  automatically  satisfied.  We  shall  address  those  conditions  in  (4.23)  that  pertain  to 
(0,b)  and  [b,°°)  after  b  has  been  determined. 

Finally,  we  note  from  (4.32)  that,  in  order  to  discuss  the  boundary  condition  (4.21)  at 
£=+oo,  V'e  must  know  the  sign  of  the  real  parts  of  r  j  and  r2 .  These  signs  depend  critically  on 

whether  the  traveling  wave  is  moving  subsonically  (1st  <  c3)  or  supersonically  (s  >  c3).  In 
continuing  the  analysis,  we  treat  these  two  cases  separately. 

5.  The  augmented  theory:  subsonic  case.  Suppose  first  that  the  given  strains  y  at 
infinity  are  such  that  the  propagation  speed  s  delivered  by  (3.6)  is  subsonic:  s  <  c3  .  Then  by 
(4.33),  (4.34),  r}  and  r2  are  both  real,  and  ^  >  0,  r2  <  0 .  From  (4.32),  it  then  follows  that  yfe) 
satisfies  the  boundary  condition  yi+°°)  =  y  if  and  only  if  the  coefficient  of  explr^-b)]  in  (4.32) 
vanishes.  Making  use  of  (4.35),  we  may  write  this  condition  as 

D,qi  +  D2q2  eq2b  =  r2(Ym  -  f) .  (5.1) 

Observe  that  both  (4.36)  and  (5.1)  must  hold  with  Dj  and  D2  given  by  (4.30).  Thus  if  one  of 

these,  say  (4.36),  serves  to  determine  b,  the  remaining  one  -  (5.1)  -  provides  a  further  necessary 

+ 

restriction  on  the  given  data  y. 

With  (5.1)  in  force,  only  the  exponentials  involving  r2  remain  in  (4.32),  so  that  y(£) 
increases  monotonically  from  y  at  §  =  b  to  f  at  £  =  +«>,  thereby  satisfying  (4.23),  and  thus 


-21- 


assuring  that  y(%)  is  indeed  in  phase  3  for  all  £  >  b. 


After  some  algebra  that  makes  use  of  (2.7),  (4.9),  (4.25),  and  (4.27)-(4.33),  one  finds  that 
(4.36)  and  (5.1)  are  equivalent  to  the  following  pair  of  equations: 


ovb/V  /Q-P\/Q-M 

Iq  +  pMq  +  r/’ 


(5.2) 


(5.3) 


where  Q,  P,  R  and  yQ  are  given  by  (4.29),  (4.25)2,  (4.34)  and  (4.27),  respectively.  We  regard 
(5.2)  as  an  equation  for  b ,  so  that  (5.3)  becomes  a  restriction  on  the  data. 

To  analyze  (5.2),  we  must  distinguish  between  real  and  complex  Q.  It  follows  from 
(4.29)  and  the  fact  that  Isl  <  c3  that  Q  is  complex  and  different  from  zero  if  and  only  if 

Islcc,*  ,  (5.4) 


where 


c**  =  c**(co)  s  { 


[co2/(l-to2)]1/2c2  <c3  if  0  <  co2  <  c2/(c2  +  c2)  , 
c3  if  to2  >  c2/(c2  +  c2) . 


(5.5) 


Figure  4  shows  the  regions  in  the  y,  f-plane  that  correspond  to  real  and  to  complex  Q. 

We  now  consider  (5.2).  Suppose  first  that  Q  is  real,  and  therefore  positive.  In  this  case, 
the  left  side  of  (5.2)  increases  with  b  from  the  value  1  at  b=0  to  +°°  at  b=+«>,  but  the  right  side  of 
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(5.2)  is  less  than  unity.  It  follows  that  (5.2)  cannot  hold  for  any  b  >  0  when  Q  is  real.  Thus  any 
point  in  either  of  the  unshaded  areas  in  Figure  4  corresponds  to  a  pair  of  strains  y,  ^  that  cannot 
be  connected  by  a  subsonic  traveling  wave  of  the  assumed  form. 

By  (4.29),  if  Q  is  complex,  it  is  pure  imaginary ,  and  both  sides  of  (5.2)  represent 
complex  numbers  on  the  unit  circle.  It  follows  that  (5.2)  has  infinitely  many  roots.  Let  0p  and 

0R  stand  for  the  respective  arguments  of  the  complex  numbers  P  +  Q  -  P  +  ilQI  and  R  +  Q  = 

R  +  ilQI.  Since  P,  R  and  IQI  are  positive,  0p  and  0R  may  both  be  taken  in  (0,  n/2).  The  positive 
roots  of  (5.2)  are  given  by 


b  =  bn=^ioi[2(n+1)7C  *2(0p  +  9r>]»  n  =  0,1,2,.... 


(5.6) 


It  can  be  shown  that,  of  all  the  roots  bn  listed  above,  only  b  =  bQ  furnishes  an  interval  (0,b)  on 
which  y(t,)  as  given  by  (4.26)-(4.31)  is  everywhere  in  phase  2;  indeed,  when  b=bQ,  y(!;)  increases 
monotonically  on  [0,b],  Thus  (5.2)  determines  exactly  one  acceptable  value  of  b. 

Having  found  b,  we  turn  to  (5.3).  Since  s  and  yQ  are  determined  by  the  given  data  y,  y 
through  (3.6),  (4.12),  (4.27)  and  (4.9),  and  since  P,  R  and  b  are  determined  in  terms  of  s  and 
material  parameters,  it  follows  that  (5.3)  represents  a  restriction  on  the  data.  In  order  to  interpret 
this  restriction,  we  first  note  that  (5.3)  can  be  solved  for  yQ  in  terms  of  s: 


(5.7) 


From  (2.7),  (4.9)  and  (4.27),  one  can  show  that 
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(5.8) 


In  (5.11),  R  and  P  are  given  by  (4.34)j  and  (4.25)2 ,  respectively,  and  b  is  understood  to  be  bQ; 
see  (5.6).  Equation  (5.9)  represents  a  restriction  on  the  given  data.  Since  by  (4.25)2,  (4.29)2  and 
(4.34)^  P,Q  and  R  depend  on  the  constants  X  and  v  of  the  augmented  theory  only  through  the 
parameter  co  given  by  (2.9),  it  follows  from  (5.6),  (5.1 1)  and  (5.10)  that  vbA,  H  and  G  have  this 
property  as  well. 

In  summary,  we  have  shown  that  if  there  is  a  subsonic  traveling  wave  of  the  assumed 
form  connecting  a  phase- 1  strain  yat  £  =  to  a  phase-3  strain  f  at  £  =  +°°,  and  if  neither  Q  nor 
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R  vanishes,  then  y(%)  must  be  given  by  (4.24),  (4.26)  and  (4.32),  b  must  be  given  by  (5.6)  with 
n=0,  and  the  given  data  y,  $  must  satisfy  the  restriction  (5.9). 

Conversely,  if  the  data  fulfill  (5.9)  and  if  b  is  given  by  (5.6)  with  n=0,  then  neither  Q  nor 
R  vanishes,  and  y(£)  as  given  by  (4.24),  (4.26)  and  (4.32)  is  a  subsonic  solution  of  the  assumed 
form  to  the  traveling  wave  problem  (4.8)-(4.10). 

+ 

6.  The  augmented  theory:  supersonic  case.  Next,  let  the  given  strains  y  at  infinity  be 
such  that  s  is  supersonic:  c3  <  s  <  c* .  In  contrast  to  the  subsonic  case,  ^  and  r2  may  now  be 

either  real  or  complex.  However,  according  to  (4.33)  and  (4.34),  when  ^  and  r2  are  real,  both 
are  negative,  and  when  they  are  complex,  both  have  negative  real  parts.  It  then  follows  that,  in 

either  case,  the  restriction  of  y(£,)  to  [b,  °°)  as  given  by  (4.32)-(4.34)  satisfies  the  boundary 

+  i 

condition  y(+°o)=y  without  further  restriction  on  the  data  y. 

We  now  return  to  (4.36)  in  order  to  determine  b  in  the  supersonic  case.  Detailed  analysis 
of  (4.36)  shows  that,  when  qj  and  q2  are  both  real,  (4.36)  has  exactly  one  positive  root  b,  and  this 

root  is  such  that  the  strain  yCH,)  associated  with  the  restriction  of  y  to  (0,  b)  is  indeed  always  in 
phase  2,  as  required.  One  can  also  show  that,  when  qj  and  q2  are  complex,  (4.36)  has  infinite’y 
many  positive  roots  b,  but  only  the  smallest  of  these  leads  to  strains  y(4)  that  belong  to  phase  2 
for  0  <  £  <  b.  Thus  whether  and  q2  are  real  or  complex,  (4.36)  determines  exactly  one 

acceptable  value  of  b  in  the  supersonic  case.  Finally,  it  can  be  shown  that,  for  this  value  of  b,  the 
strain  y(^)  associated  with  the  restriction  (4.32)  of  y  to  [b,  °°)  is  in  phase  3  for  all  %  >  b  in  those 
supersonic  cases  for  which  ^  and  r2  in  (4.33)  are  real.  We  have  been  unable  to  show  that  this 
conclusion  persists  in  supersonic  motions  with  complex  r’s,  but  numerical  calculations  suggest 
that  it  does.  We  shall  take  this  for  granted  in  what  follows. 
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In  summary,  if  there  is  a  supersonic  traveling  wave  of  the  assumed  form  connecting  a 
phase- 1  strain  y  at  £,  =  -«>  to  a  phase-3  strain  ^  at  ^  =  +°°,  and  if  neither  Q  nor  R  vanishes,  then  a 
unique  value  of  b  can  be  found,  and  y(£)  must  be  given  by  (4.24),  (4.26)  and  (4.32).  Conversely, 
given  ye  (-1,0)  and  f  >  Ym  (see  Figure  2)  for  which  Q  and  R  do  not  vanish,  then  for  a  suitable 

value  of  b,  y(^)  as  given  by  (4.24),  (4.26)  and  (4.32)  represesents  a  supersonic  traveling  wave 

+ 

without  further  restriction  on  the  data  y. 

Finally,  we  note  that  the  special  case  of  the  augmented  theory  for  which  X=0,  v  >  0 
corresponds  to  a  simple  viscoelastic  model  of  the  bar.  Although  we  do  not  treat  this  case  in  detail 
here,  we  remark  that  it  leads  to  the  conclusion  that  all  traveling  waves  in  the  augmented  theory 
are  supersonic.  A  related  observation  was  made  by  Slemrod  [5]  in  motivating  the  inclusion  of 
capillarity  as  well  as  viscosity  when  augmenting  the  theory  of  the  van  der  Waals  fluid.  As  noted 
in  [4],  [17],  slowly  moving  phase  boundaries  in  solids  have  been  observed  in  experiments;  the 
special  case  X=0  of  the  present  augmented  theory  would  therefore  seem  too  special. 

At  the  other  extreme,  the  special  case  of  the  augmented  theory  in  which  v=0,  X>0  leads  to 
purely  dissipation-free  behavior,  as  suggested  by  (4.11). 

7.  Permissible  phase  boundaries  ;n  the  elastic  bar.  According  to  the  criterion  based 
on  the  augmented  theory  described  here,  t  sible  phase  boundaries  within  the  elastic  theory 
are  those  for  which  the  associated  strains  y  and  f  can  be  connected  by  a  traveling  wave  in  the 
augmented  theory.  For  phase- 1,  phase-3  interfaces,  pairs  of  strains  y,  y  permitted  by  the  jump 
conditions  (2.4),  (2.5)  in  the  elastic  theory  correspond  to  points  in  the  region  F  of  Figure  2. 
Among  all  points  in  T,  those  for  which  -1  <  y  <  0  correspond  to  supersonic  phase  boundaries, 
and  those  with  0  <  y  <  yM  are  associated  with  subsonic  phase  boundaries. 
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As  we  have  shown  in  the  preceding  section,  every  "supersonic  pair"  y,  (except  possibly 
those  for  which  either  R  or  Q  vanishes)  can  be  connected  by  a  traveling  wave  in  the  augmented 
theory.  Thus  all  such  supersonic  pairs  are  "preferred"  according  to  the  present  criterion. 

In  contrast,  the  results  of  Section  5  show  that,  among  all  points  (y,  y)  in  T  that  correspond 
to  subsonically  propagating  phase  boundaries,  only  those  that  satisfy  the  restriction  (5.9)  are  such 
that  yean  be  connected  to  y  by  a  traveling  wave  in  the  augmented  theory.  Interpreted 
geometrically,  condition  (5.9)  is  the  parametric  representation  of  a  curve  in  the  subsonic  portion 
of  T  that  corresponds  to  the  locus  of  pairs  y,  that  are  preferred  in  the  elastic  theory  according  to 
the  present  criterion. 

More  significantly,  the  restriction  (5.9)  on  preferred  .ubsonic  phase  boundary  strains  is 
equivalent  to  a  kinetic  relation  of  precisely  the  form  (3.12).  By  substituting  for  y  and  y  from 
(5.9)  into  the  representation  (3.9)  for  the  driving  traction  f  at  a  phase- 1,  phase-3  interface,  we 
find  indeed  that 


f  =  <p(s)  , 


where 


(7.1) 


cp(s)  =  -  4ij  -  H3) 


G2(s) 


,2  -2w  2 

(Cj  -  S  )(Cj  ■ 


S2) 


]■ 


"C**  <  s  <  c,* 


(7.2) 


In  (7.2),  G(s)  is  given  by  (5.10),  (5.11). 

It  is  appropriate  to  note  here  that  the  solution  y(%)  of  the  traveling  wave  problem  in  the 
augmented  theory  has  a  limit  as  the  viscosity  v  tends  to  zero,  the  parameter  £0  and  the  data  y,  y 
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being  held  fixed.  Because  P,  Q  and  R  depend  only  on  co  and  the  data,  it  follows  from  (5.6)  that 
the  length  b  of  the  phase-2  interval  tends  to  zero  in  this  limit.  It  is  also  readily  shown  that  y(£) 
tends  to  y  for  every  £<0  and  to  ^  for  every  ^>0  as  v  ->  0,  co  fixed.  This  zero-viscosity  limit  of 
y(£,)satisfies  all  the  differential  equations  and  jump  conditions  of  the  elastic  theory.  Within  that 
theory,  it  represents  a  propagating  phase  boundary  moving  with  a  velocity  s  that  satisfies  (3.8) 
and  the  entropy  admissibility  requirement  (3.5).  Moreover,  because  G  in  (5.10)  depends  on  v  and 
X  only  through  0),  the  kinetic  relation  (7.2)  remains  unaffected  by  the  limit  process. 

Figure  5  describes  the  special  kinetic  response  function  cp  of  (7.2)  by  plotting  f/fQ  vs.  s/c3 
according  to  (7.1),  (7.2);  here  fQ  =  (|ij  -  |i3)yMYm/2 .  One  finds  that  f/fQ  depends  only  on  the 
three  material  parameters  oM/om,  y  /y  and  CD;  the  figure  is  plotted  for  fixed  values  of  the  first 
two  of  these  and  for  various  values  of  to.  For  small  values  of  co  (corresponding,  for  example,  to 
large  viscosity),  the  graph  strikes  the  top  and  bottom  boundaries  of  the  permissible  region  in  the 
s,  f-  plane  when  s=c*„  .  This  occurs  when  co  <  c3/(c2  +  c3)  ;  see  (5.5).  For  the  other  values  of 

co  considered,  f/fQ  approaches  -«>  as  s/c3  tends  to  -1,  but  f/fQ  tends  to  unity  as  the  latter  ratio 

approaches  +1.  For  large  values  of  co,  or  small  viscosity,  the  kinetic  relation  corresponds  to 
nearly  dissipation-free  motion  of  the  phase  boundary  over  most  of  the  range  of  velocity  s.  Note 
that,  from  the  graphs,  cp(s)  increases  monotonically  with  s  for  the  values  of  the  parameters 
considered.  The  kinetic  relations  encompassed  by  the  analysis  in  [1]  were  assumed  to  have  this 
monotonicity  property.  It  can  be  shown  that,  as  indicated  in  Figure  5  and  as  expected  from  (4.1 1) 
and  (4.12),  <p(s)  as  given  by  (7.2)  satisfies  the  requirement  (3.13)  imposed  on  every  kinetic 
response  function  by  the  entropy  admissibility  condition  (3.5). 
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correspond  to  a  special  case  of  those  that  have  the  form  (3.12)  and  were  considered  in  [1,3]. 
These  in  turn  are  by  no  means  the  most  general  kinetic  relations  conceivable.  Moreover, 
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augmentation  is  not  the  only  way  to  specify  the  kinetics  of  a  phase  transition  within  a  continuum 
model;  thermal  activation  theory,  for  example,  provides  an  entirely  different  approach  that  may 
be  ch}  sic  .illy  appropriate  under  certain  circumstances.  Even  within  the  general  framework  of 
auc  n  miation,  alternatives  suggest  themselves:  one  might  supplement  the  Newtonian  viscous 
dai?  ■  g  vyt  in  (2.8)  by  adding  a  term  proportional  to  y^;  the  resulting  detailed  kinetics  would 

presumably  differ  from  those  described  by  (7.1)  and  (7.2).  Indeed,  it  is  unlikely  that  any  scheme 
base  on  an  augmented  theory  for  modeling  the  kinetics  of  phase  transitions  in  solids  can  lay 
claim  io  universality. 
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FIGURE  1.  STRESS-STRAIN  CURVE  FOR  THE  TRILINEAR  MATERIAL. 


FIGURE  2.  TOE  REGION  T  IN  TOE  y ,  y  -  PLANE. 


FIGURE  3.  ADMISSIBLE  IMAGE  OF  T  (INSET)  IN  s  ,  f  -  PLANE . 
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FIGURE  4.  SUBREGIONS  OF  T  CORRESPONDING  TO  REAL  AND  COMPLEX  Q . 
L:  y  =  Ay,  A=  [(1-  cd2)  c^-a2^]  /  [(1  -o>2)  c|  -a?c}] . 

Q  IS  COMPLEX  IN  HATCHED  REGIONS. 
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